*******  replication Table A12: negative binomial models

eststo clear


** MNLF
use "MNLF_Target_Final.dta" , clear


quietly: nbreg attack_region ramadan christian_4days christian_2days christian_1days nonreli_holiday  /*
*/ nation_election Islamic_3day Islamic_2day Islamic_1day  rain hot lag_attack_region cubic1-cubic3 i.region  , vce(cluster region) from (ramadan=-1.1)

eststo model1


** MILF: clogit does not work; ombinations results in numeric overflow; computations cannot proceed

use "MILF_Target_Final.dta" , clear

quietly: nbreg attack_region ramadan christian_4days christian_2days christian_1days nonreli_holiday  /*
*/ nation_election Islamic_3day Islamic_2day Islamic_1day  rain hot lag_attack_region lag_gov_counter cubic1-cubic3 i.region  ,  vce(cluster region) from (ramadan=-0.5)
eststo model2



** BIFF
use "BIFF_Target_Final.dta" , clear

quietly: nbreg attack_region ramadan christian_4days christian_2days christian_1days nonreli_holiday  /*
*/ nation_election Islamic_3day Islamic_2day Islamic_1day  rain hot lag_attack_region lag_gov_counter cubic1-cubic3  i.region  , vce(cluster region) from (ramadan=0.65)
eststo model3



*** Patani before 2000
use "Patani_bef_Target_Final.dta" , clear

quietly: nbreg attack_region ramadan  Bud_holiday nonreli_holiday  /*
*/ nation_election Islamic_3day Islamic_2day Islamic_1day rain hot  lag_attack_region cubic1-cubic3  lag_gov_counter  i.region,    vce(cluster region) from (ramadan=-0.05)
eststo model4


*** Patani after 2000
use "Patani_aft_Target_Final.dta" , clear

quietly: nbreg attack_region ramadan  Bud_holiday nonreli_holiday  /*
*/ nation_election Islamic_3day Islamic_2day Islamic_1day rain hot  lag_attack_region cubic1-cubic3 lag_gov_counter   i.region, vce(cluster region) 

eststo model5



*** GAM
use "GAM_Target_Final.dta" , clear

quietly: nbreg attack_region ramadan christian_4days christian_2days christian_1days nonreli_holiday  /*
*/ nation_election Islamic_3day Islamic_2day Islamic_1day  rain hot lag_gov_counter lag_attack_region  cubic1-cubic3 i.provincode, vce(cluster provincode)  from (ramadan=-0.1)
eststo model6


esttab using Table_A12.tex,  order(ramadan Islamic_3day Islamic_2day Islamic_1day  christian_4days christian_2days christian_1days Bud_holiday    nonreli_holiday  nation_election lag_gov_counter  rain hot attack_gov_region_lag cubic1 cubic2 cubic3 ) label star(+ 0.1 * 0.05 ** 0.01 *** 0.001) se b(3) replace booktabs ///
  alignment(D{.}{.}{-1})      mtitles("Model 1" "Model 2"  "Model 3" "Model 4"  "Model 5" "Model 6")                   ///
  title(Fixed Effects Model table\label{tab1})
    
  
*******  replication Figure A2: negative binomial models
  
coefplot model1 model2  model6 model4 model3   model5 , keep(ramadan ) xline(0) order(ramadan ) saving(Figure_A2)

